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SUMMARY 

In  this  paper  we  consider  the  time-dependent  behaviour  of  a single-server 
queue  with  Poisson  input  and  constant  service  time,  where  the  input  rate  equals 
the  reciprocal  of  the  service  time.  Formulae  for  the  probabilities  associated 
with  different  queue-lengths  have  been  derived  by  an  empirical  method  which 
included  inspection  of  some  numerical  results.  These  formulae  are  now  proved 
correct  by  induction. 
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I INTRODUCTION 

The  problems  studied  in  this  paper  were  raised  by  Attwooll  when  preparing 
a lecture*  on  time-dependent  queues.  It  is  well  known  that  a single-server  queue 
with  Poisson  arrivals  and  constant  service  time  (an  M/D/ I queue)  lias  a steady- 
state  queue-length  distribution  if  and  only  if  the  arrival  rate  is  less  than  the 
service  rate.  (The  service  rate  is  the  reciprocal  of  the  service  time.)  If  the 
arrival  rate  equals  (or  is  greater  than)  the  service  rate  the  queue  is  'unstable' 

and  the  average  queue  length  increases  indefinitely  with  time.  The  steady-state 
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solution  for  the  M/D/ 1 queue  was  published^  in  1909.  The  time-dependent  solution 
for  the  M/M/ 1 queue  (exponential  distribution  of  service  times)  was  first  published 
in  1953  3. 

In  this  paper  we  consider  the  case  where  the  arrival  rate  is  exactly  equal 
to  the  service  rate.  Without  loss  of  generality  we  define  unit  time  to  equal  the 
service  time  (therefore  equal  also  to  the  mean  time  between  arrivals). 

Let  q(n)  be  the  probability  that  there  are  i customers  in  the  queue  at 
time  n . We  shall  only  consider  integer  values  of  n . The  "number  in  the 
system"  includes  the  customer  being  served,  if  there  is  one. 

The  number  of  arrivals  during  a unit  interval  has  a Poisson  distribution 
with  mean  I,  so  the  probabilities  of  0,  1,  2,  j arrivals  are  respectively 
exp(-l),  exp(-l),  exp(-!)/2!,  exp(-l)/j!  . The  number  of  departures  during  a 
unit  interval  is  I (if  there  is  a customer  in  the  system  at  the  start  of  the 
interval)  or  0 (if  the  system  is  empty  at  the  start  of  the  interval).  It  is 
easy  to  write  down  equations  giving  the  queue-length  probabilities  at  time  (n  + 1) 
in  terms  of  the  probabilities  at  time  n : 


Vn  + 

1) 

- |qQ(n)  + q j (n) } exp(-l) 

(1) 

q j (n  + 

1) 

= |qQ(n)  + q j (n)  + q2(n)}  exp  (—  1 ) 

(2) 

<l2(n  ♦ 

1) 

- Jq^(n)/2!  + q)(n)/2!  + q2(n)  + q.j(n)|  exp(-l) 

(3) 

q.(n  + 1)  - ^qQ (n) /^^ ! + q ] (n) / £ ! + q2(n)/(i  - I)!  + ...q^in)  + qi+,(n)|  exp(-l). 

(4) 

We  consider  the  case  where  the  system  starts  empty,  so  q^(0)  = I and  all 
the  other  initial  probabilities  equal  0 . Values  of  the  different  probabilities 
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for  n = 1,  2,  3,  etc  can  then  be  obtained  numerically  by  repeated  solution  of  the 
above  set  of  equations,  using  a computer.  Attwooll  did  this  and  noticed  that  the 
probability  of  no  customers  waiting  was  given  by  the  formula: 

( + nn+1 

q0(n)  + q((n)  = ~~('n~T  1) ! exP  (~  n)  • (5) 

He  was  unable  to  prove  this  formula.  We  shall  proceed  by  assuming  its  correctness 
in  order  to  deduce  results  for  other  queue-lengths,  then  verify  our  results  by 
induction. 

2 DERIVATION  OF  FORMULAE 

Using  (5)  to  substitute  in  equation  (1): 

( + n n+ 1 

q (n  +1)  = — — exp  (-  n - 1) 

(n  + 1)! 


so  obviously 


q0(")  = exp(-  n) 


(6) 


This  expression  is  undefined  for  n = 0.  We  assume  q^CO)  “I  , ie  the  system 
starts  empty. 

To  simplify  the  equations  we  introduce  the  notations 


Qi(n)  = q. (n)  exp(n)  (7) 


and 


Ini  = nn/n! 


|0)  = 1 


so  equation  (6)  can  be  put  in  the  form 


Qg(n)  = Ini 


(8) 


From  equations  (1)  to  (4)  we  can  derive  the  system  of  equations: 
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qQ(n  + I) 

Qj (n  + I) 

Q,(n  + 1) 


Vn)  * ‘'i 

(n) 

(9) 

V">  * <*1 

(n)  + Q2(n) 

(10) 

Q0(n>/2!  . 

q,(n)/2!  + Q2(n)  + q^n) 

(II) 

Qj(n  + I)  = Q0(n)/j!  + Q , (n) / j ! + Q2(n)/(j  - 1)!  + ...  + (n)  + Q(^+1)(n). 

(12) 

Using  equation  (8)  to  substitute  for  both  QQ(n)  and  g^Cn  +1)  in 

equation  (9),  and  re-arranging  terms,  we  find: 

Q,(n)  = [n  + 1]  - In)  (13) 

Using  equations  (8)  and  (13)  to  substitute  for  Q((n  + 1)  , Qj(n)  and 
Q^(n)  in  equation  (10): 


Q2(n)  = (n  + 2]  - 2 [n  + l]  (14) 

Proceeding,  in  the  same  way,  to  substitute  into  equation  (11),  we  find: 


Q-j(n) 

= 

In  + 

3| 

- 3 In 

+ 2]  + | 

In  + 

1] 

(15) 

and  it 

can  also  be 

shown  that: 

Vn) 

- In  + 4| 

- 4 In  + 

3]  + 

4.2 

2! 

In 

+ 2]  - 

jr  |n  + 

1) 

(16) 

Q5(n) 

“ In  + 5l 

- 5 In  + 

4)  ♦ 

5.3 

2! 

In 

+ 3)  - 

5.22  , 

3!  In 

+ 2]  ■ 

In  + 

II 

(17) 

Vn) 

• In  ♦ 6| 

- 6 In  + 

5)  + 

6.4  In  •» 
2! 

• 41  _ 

6.32  . 

3!  ln  + 

3]  + 

6.23 

4! 

(n  + 

2]  • 

6 In  ♦ 1 1 
5! 

• 

• 

• 

(18) 

g.  (n) 

i-l 

' ‘‘I 

(j  - i)j 

j! 

-1 

— In 

i ♦ i 

- 

jl 

(i  > 

1)  . 

(19) 

j-0 
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The  distribution  of  queue  lengths  can  now  be  calculated,  using  equation  (7)  to 
transform  equations  (8),  (13)  and  (19)  back  into  the  original  variables  q.  : 

n 

q (n)  = exp(-  n)  = Ini  exp(-  n) 

q,<n>  = V+'jyr  " ^rj  exp(~ n)  = (l>  + ^ " M)  exp(" n)  (: 


q^(n)  = - i exp (-  n)  ^ 

3=0 


(j  ~ i)^"1  tp  ♦ i ~ i 


(i  > 1)  • 


By  assuming  the  correctness  of  Attwooll's  empirical  formula  (equation  (5)) 
we  have  generalised  his  result  and  obtained  formulae  (equations  (20),  (21),  (22)) 
for  all  queue  lengths. 

3 PROOF  OF  CORRECTNESS 

We  shall  prove  these  formulae  by  induction,  ie  we  shall  prove  first  that  the 
formulae  are  correct  for  n = 0,  and  secondly,  provided  the  formulae  are  true  for 
all  times  up  to  and  including  some  value  n they  are  also  correct  for  time 
(n  ♦ 1). 

For  n = 0,  straightforward  substitution  in  equations  (20)  and  (21)  gives: 


qQ (o)  = i , q , (0)  = o 


Putting  n = 0 in  equation  (22)  gives: 


q j (o) 


i V u - « - J)1' 

1 L J!  a - j>! 


(i  > 1) 


hence 


V0)  • a-hyr  £ <-  l,J<i " j)*''  (i) 


writing  k for  (i  - j)  and  summing  in  the  reverse  order: 
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(i  > 1)  . 


qi(0) 


k=  I 


1)  k 


In  the  Appendix  it  is  shown  that 


I 


(-  0 


k,  i-  I 


0 


so  that 


q • (0)  = 0 

i 


(i  > 1)  . (24) 


Together  with  equations  (23)  this  completes  the  first  part  of  the  proof. 

Now  assuming  equations  (20)  to  (22)  are  correct  for  some  value  n , the 
probabilities  at  time  (n  + 1)  are  found  by  substituting  for  the  q^(n)  on  the 
right-hand  side  of  equations  (1)  to  (4).  Thus  we  find  that 


qQ(n  + I) 


(n  + 1) 


n+ 1 


(n  + D! 


exp(-  n - 


1) 


and  this  is  obviously  of  the  correct  form, 
q((n  + I)  = (|nl  + In  + 1 1 - [n]  + In 
= |[n  + 2]  - [n  + l])  exp (- 


so  equation  (20)  holds  at  time  (n  + I). 

+ 2]  - 2 In  + l]j  exp(-  n - 1) 
n - I ) 


and  this  is  also  of  the  correct  form,  so  equation  (21)  holds  at  time  (n  + 1). 

Equation  (3)  is  just  a special  case  of  equation  (4)  which  can  be  written  as 


q. (n  + 1) 


q (n)  q (n) 

-~r-. + — : 


i+1 


. V Vn> 

/ ( i 4.  1 - 1 U 


exp (-  1) 


(25) 


Using  equations  (20)  to  (22)  to  substitute  for  the  q (n)  on  the  right-hand 
side  of  equation  (25) 


q^(n  + I)  = exp (-  n - 1) 


i+ I i-1 


(nl  [n  + ll  ~ [nj 

tit;  <u! 


.vk-l 

j In  + j - kl 


j=2  k 


q. (n  + 1)  - exp  (-  n - 1 ) 


i+1  j=l 


ZZ^ 

j-1  k=0 


.,k-l 
. ~ J) 

J k! 


In  + j - k] 


(26) 


Now  let  m=i+l-j+k  so  that  j=i+k+!-m.  To  transform  the 
summation  it  is  helpful  to  draw  up  a table,  showing  which  combinations  of  j and 
k are  to  be  summed  over,  and  the  value  of  m for  each  summable  combination 
see  Fig  1.  Substituting  for  j , and  summing  over  m and  k instead  of  over 
j and  k : 


q^(n  + 1) 


= exp(-  n - 1) 


l m 

r-'  v- » 


zz 


(i  + k + 1 - m)  (m  - i - 1 ) 


k-I 


k! 


In  + I + i - m] 


m=0  k=0 


q^(n  +1)  = exp(-  n - 1) 


u 


m . ^ 

(i  + k + I - m)  (m  - i - I ) f 

i. -n  * 


In  + 1 + i - m] 
m-0  ' k=0 


(27) 


but 


m 


z 


(i  + k+l-m)(m-i-l) 


k-1 


m 


k! 


-z<- 


I) 


k-1  (i  + 1 - m + k)  (i  + I - m) 


k-1 


k! 


k=*0 


k-0 


Z <- ,,k" 


-I  I (i  + 1 - m)k  A (i  + I -m^Vk^l 


k! 


(k  - 1 ) 


k-0 


(-  1) 


m- 1 (i  + 1 - m) 


m 


ml 


(28) 


substituting  this  last  expression  for  the  inner  sum  in  equation  (27)  we  obtain: 
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q^(n  ♦ 1)  = 


exp (-  n - 1 ) ) (-  1 ) 

m=0 


1 

2> 


m (m  - i - l)m  * [n  + I + i - m] 
m! 


(29) 


By  comparing  equations  (29)  and  (22),  it  can  be  seen  that  the  formula  is  correct 
for  time  (n  + 1).  This  completes  the  second  part  of  the  proof. 

We  have  now  shown  that,  provided  equations  (20),  (21)  and  (22)  are  correct 
up  to  time  n , they  are  also  true  for  time  (n  + 1).  But  we  showed  earlier  that 
the  equations  are  correct  for  time  0 , therefore  they  are  also  correct  for  time  1, 
therefore  they  are  also  correct  at  time  2,  and  so  on.  Thus  equations  (2)  to  (22) 
are  true  for  all  integer  values  of  the  time  n . 


We  have  thus  proved  the  correctness  of  equations  (20)  to  (22).  Attwooll's 
empirical  result,  equation  (5),  can  easily  be  proved  by  adding  equation  (20)  to 
equation  (21). 

4 MEAN  QUEUE  LENGTH 

The  mean  queue  length  may  be  defined  as 


q (n)  = ^ (i  ” l)qi(n)  (30) 

i=2 

(the  customer  being  served  is  not  included).  By  multiplying  each  of  equations  (1) 
to  (4)  by  an  appropriate  factor,  and  summing  over  i , Attwooll  obtained  the 
result 

q (n  + 1)  = q (n)  + {qQ(n)  + qj(n)}e  1 . (31) 

Using  equation  (5),  which  we  have  proved  in  section  3 

(32) 

(33) 

For  large  values  of  m , the  value  of  m!  is  given  approximately  by  Stirling's 
formula 


vn+l 


q(n  + 1)  - q (n) 


(n  + 1 ) z , s 

Tn'~)T  6XP(~  n ' 0 


hence 


Em 

exp  (-  ra) 

m=  1 


■ 
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m!  * (27rm)^‘^Dim  exp(  - m) 


Thus  for  large  values  of  n 


q(n)  - q (n  - 1)  = (2irn) 


-0.5 


Hence  it  may  be  shown  that 


q( 


-> 


(34) 


(35) 


Comparing  approximate  numerical  values  obtained  from  equation  (35)  with  the 
corresponding  exact  values  obtained  from  equation  (33),  it  was  found  that  a very 
good  approximation  is  given  by  the  equation 


q(n) 


2n  + 

IT 


1.2 


0.66  . 


(36) 


5 SUMMARY  OF  RESULTS 

The  following  results  (except  for  the  last  equation  (A-2))  relate  to  the 
queue  system  M/D/1  with  arrival  rate  equal  to  service  rate. 

The  probability  of  i customers  in  the  system  (including  the  customer 
being  served)  at  time  n is  given  by  the  formulae 


q0(n)  = 

[n]  exp(-  n) 

(20) 

q,(n)  = 

( [n  + 1]  - [n])  exp  (-  n) 

(21) 

q{(n)  = 

- i exp (-  n)  £ (j-~ 

(22) 

i=o 

where  [rl  » rr/r! 
and  [0]  - 1 . 

The  expected  queue-length  (excluding  the  customer  being  served)  at  time  n 
is  given  by  the  formula 
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Appendix 

SUMMING  A FINITE  SERIES 


In  this  Appendix  we  find  the  sum  of  the  series: 


S(n,r) 


z 


(-  l)kk 


(A— 1 ) it  follows  that 


OO 


Z 

r-0 


S(n,r)zr 

r! 


By  equating  coefficients  of 


r 


z 


it  can  be  seen  that 


S(n,r)  =0  0 $ r < n 


y <-  o\r(?\ . 


0 


0 % r < n 


The  special  case  S(n,l)  = 0 is  included  in  Ref  4,  on  page  4 as  formula 
0. 154.2. 


(A-3) 


(A-4) 


(A-5) 
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